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Intermediate state (IS) flux structures in mesoscopic type-I superconductors are investigated 
within the Ginzburg-Landau theory In addition to well-established tubular and laminar struc- 
tures, the strong confinement leads to the formation of (i) a phase of singly quantized vortices, 
which is typical for type-II superconductors and (ii) a ring of a normal domain at equilibrium. The 
stability region and the formation process of these IS flux structures are strongly influenced by the 
geometry of the sample. 
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Intermediate state (IS) of type-I superconductors has 
received a revival of interest in recent years [TJ [2] as one of 
the rare systems where competing interactions lead to the 
formation of spatially modulated structures. This kind 
of complex structures have been observed in other sys- 
tems like ferrofluids, amphiphilic monolayers and chem- 
ical reaction-diffusion systems [3]. Analogies to type-I 
superconductors also extend to astrophysics [4]. 

IS flux structures can be described theoretically as- 
suming parallel stripes of normal and superconducting 
(SC) domains [5] or a periodic array of multiquanta flux 
tubes [B]. The tubular structure is more mobile, while the 
laminar structure is topologically constrained [13 and 
transition between these states occurs for larger volume 
fraction of normal domains [14]. However, these regu- 
lar patterns are rarely seen in experiment [7 . Instead, 
much more complex flux structures - highly branched 
and intricate fingered patterns of flux domains - are often 
observed El GB [IS EE] • In addition, type-I super- 
conductors show clear hysteresis - flux tubes are obtained 
during magnetic field penetration and laminar structures 
are formed during magnetic field expulsion [TTJ [12] , sug- 
gesting that the system gets stuck in a metastable state, 
i.e., the sample is not in thermodynamic equilibrium. 
Zero-bulk pinning [2] superconducting spheres and cones 
show no hysteresis with flux tubes dominating the IS. 
With increasing applied magnetic field, a new phase - a 
suprafroth - is formed [lj, which represents neither tubu- 
lar nor laminar structure. 

The behavior becomes richer in the mesoscopic regime, 
where, in addition to the competition between the mag- 
netic energy that favors the formation of small normal 
domains and the positive surface energy that tends to 
form large domains, confinement effects become impor- 
tant. In type-II superconductors confinement strongly 
influences the distribution of vortices leading to, e.g., the 
formation of giant vortices [15] and symmetry induced 
vortex-antivortex pairs [IB] . Misko et al., stabilized such 
vortex-antivortex patterns in a long type-I superconduct- 
ing prism with triangular cross section [17], due to the 



competition between vortex-antivortex repulsion and fi- 
nite size effects. However, their study was limited to a 
very narrow range of applied magnetic field. Our main 
goal in this paper is a systematic study of ground state 
flux structures in mesoscopic type-I superconductors. We 
also address another important problem not yet stud- 
ied neither experimentally nor theoretically: the effect of 
sample geometry on these mesoscopic IS patterns. 

From the theoretical point of view, the difficulty of 
modeling the IS comes from the 3D nature of the mag- 
netic interaction of the domains. Therefore, approximate 
expressions are usually used for the magnetic energy. 
[BJ [14j [18j [19] In our approach we used the Ginzburg- 
Landau (GL) theory, where no approximation for the 
magnetic energy is used and the structure of the domains 
is not predetermined. We solve the GL equations for the 
order parameter \I/ and the vector potential A: 

HV- A) 2 ^ = (1- |^| 2 )^, (1) 
-k; 2 V x V x A = Im(^*W) - |^| 2 A, (2) 

where n is the GL parameter. Here, the distance is mea- 
sured in units of the coherence length £, the vector poten- 
tial A in ch/2e£, and the order parameter \I> in ^—a//3 
with a, (3 being the GL coefficients. Following the nu- 
merical approach of Ref. [20] , we apply a finite- difference 
representation for the order parameter and the vector 
potential on a uniform 3D Cartesian space grid (up to 
256 x 256 x 256 grid points). \I/ satisfies the boundary 
condition (— iV — A)\I/| n = at the sample surface and 
far away from the superconductor A is determined by 
the external applied field. The simulations of the IS are 
conducted in i) field sweep up: we started from the full 
Meissner state (|\I/| = 1) and slowly increased the mag- 
netic field, after reaching the stationary state; ii) field 
sweep down: we started simulations with \^\ =0 and 
H > H c and decreased the field with small steps; and 
iii) field cooled simulations starting from random initial 
conditions for each value of the applied field. 

Equilibrium flux structures. - As a representative 
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FIG. 1: (Color online) The ground state flux structures for a 
rectangular prism of size a and height b (see the inset) as a 
function of a and the applied magnetic field H for b = 4£ and 
k, — 0.4. The insets (2-6) denote the Cooper-pair density plots 
[red (dark gray)/yellow (light gray) corresponds to high/low 
\ip\ 2 ] of the corresponding states. 

example of mesoscopic type-I superconductors, we stud- 
ied IS flux structures in a rectangular prism of size a and 
height b exposed to a homogeneous field H (see insets 
of Fig. [I]). We constructed the equilibrium phase dia- 
gram as a function of a and H for n = 0.4 and 6 = 4£ 
as shown in Fig. [T] The ground state was found from 
the field cooled simulations by comparing the energies of 
all found states for given magnetic field and sample pa- 
rameters. The Meissner state (region 1) is found at small 
fields and the stability of this phase increases to larger 
fields with decreasing sample size. The most prominent 
feature of this phase diagram is the existence of region 
2 where singly quantized fluxes (i.e., individual vortices) 
are nucleated in the sample (see inset 2). Note that this 
state is typical for type-II superconductors and was never 
found in bulk type-I samples [2j [TTJ [24] , where the energy 
of singly quantized vortices is larger than the energy of 
multiply quantized flux tubes. At higher fields individ- 
ual vortices are too close to each other and giant vortices 
nucleate together with individual vortices (region 3) as 
illustrated by inset 3 of Fig. [I] In region 4 laminar-like 
structures start to form in combination with flux tubes 
(inset 4). Such laminar-like structures are obtained for 
larger a, usually in the form of a ring of normal domain 
(inset 5). This novel state is a consequence of strong 
confinement imposed by the sample surface and resem- 
bles recently found suprafroth in macroscopic SC spheres 
PQ . Surface superconductivity (region 6) is formed before 
the system transits to the normal state (inset 6). Notice 
also that, we did not find branching of normal domains, 
which occurs in very thick samples [7 . From this phase 
diagram we notice that the ground state flux structures 



in mesoscopic type-I superconductors do not depend only 
on the applied magnetic field (as in the case of bulk su- 
perconductors [H]), but also on the confinement due to 
the finite size of the sample. The latter leads to the for- 
mation of singly quantized flux tubes and a ring of normal 
domain. 

Effect of sample geometry. An important fac- 
tor that influences the formation of IS patterns is the 
sample shape and geometry [2j [2TJ [26]. For example, 
a recent experiment showed [2] that zero-bulk pinning 
discs and slabs show hysteretic behavior - flux tubes ap- 
pear on magnetic field penetration and lamellae on flux 
exit. Spheres and cones show no hysteresis with flux 
tubes dominating the intermediate field region. In what 
follows, we study the influence of the sample geometry 
on the IS flux structures, as an example of a supercon- 
ducting cube and sphere with the same superconducting 
volume V. Figure [5] shows the free energy (in units of 
Fq = H^V/Stt) (a) and the magnetization (b) of a cube 
with size a — 20£ as a function of the applied magnetic 
field. In the field sweep up regime (thick red curve), the 
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FIG. 2: (Color online) The free energy (a) and the magne- 
tization (b) as a function of the applied magnetic field for a 
superconducting cube of size a = 20£ and for k — 0.4. Thick 
red (blue thin) curves are the results obtained when increas- 
ing (decreasing) the magnetic field and the symbols (hyphens) 
show the results of field cooling simulations. The inset shows 
the total number of flux in the sample as a function of the 
applied field obtained during field sweep up (filled red circles) , 
field sweep down (open blue circles) and of the lowest energy 
state from field cooled (stars) simulations. 
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sample is in the Meissner state up to H p = 0.3SH C - 
the penetration of flux is prohibited by the surface bar- 
rier. When the penetration field H p is reached four flux 
quanta enter the sample forming one circular flux tube or 
a giant vortex. With a small increase of the field, 4 flux 
quanta enter forming 4 giant vortices, each with vorticity 
2 (Fig. [3^ a)). The strong confinement in the mesoscopic 
regime prevents the formation of hexagonal structures 
and we usually obtain square symmetric structures (Figs. 
[3ja,b,e-h)). Previous studies showed that with increas- 
ing field magnetic flux enters the sample in the form of 
tubes [11] which were believed to be due to the presence 
of a surface barrier [21] . The smallness of our sample en- 
hances the effect and only a single flux quantum enters 
through each side of the cube and join the existing nor- 
mal domains. In this way, the square symmetry of the 
flux structure remains unaltered and only the number of 
flux quanta in each normal domain increases (Fig. |3jb)). 
With further increasing the field a ring of normal phase 
is formed surrounding a SC domain (Fig. [3^c)). This 
ring structure closes at H = 0.66H C and only surface 
superconductivity (Fig. [3^d)) remains until the system 
transits to the normal state. 

The topological hysteresis in the IS flux structures in 
macroscopic samples [2] was previously explained by the 
presence of a surface barrier which is absent in decreas- 
ing fields [22 - a large number of lamellae are connected 
to the sample edge during flux exit, allowing the mag- 
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FIG. 3: (Color online) Contour plots of the Cooper pair den- 
sity 2 in the middle plane of the cube obtained during mag- 
netic field sweep up (a-d), down (e-h) and in field cooling 
regime (the ground state) (i-1) at the applied fields H/H c =: 
0.44 (a), 0.52 (b), 0.55 (c), 0.66 (d), 0.43 (e), 0.37 (f), 0.15 
(g), 0.08 (h), 0.17 (i), 0.28 (j), 0.48 (k), and 0.57 (1). The 
number of fluxoids N is: 8 (a), 20 (b), 28 (c), 36 (d), 48 (e), 
40 (f), 16 (g), 8 (h), 3 (i), 8 (j), 22 (k), and 40 (1). 



netic flux to exit continuously, while these large normal 
domains are broken into smaller pieces during flux pene- 
tration. Thin blue curve in Fig. |2|a) shows the free en- 
ergy of the cube when decreasing the magnetic field H. 
Surface superconductivity first nucleates resulting in a 
narrow flux-free zone, as obtained in the experiment [23 . 
With further decreasing H, a SC phase nucleates in the 
central part of the sample, embedded within the normal 
phase (similar to the state in Fig. [3^c)) reducing the to- 
tal number of flux quanta (see the inset of Fig. [2|a)). At 
H = 0A3H C the long normal domain breaks down into 4 
smaller domains that are arranged into a square symmet- 
ric structure (Fig. [3|e)). When reducing the magnetic 
field further the length of the normal lamellae decreases 
until they are reduced to flux tubes (Fig. [3^f)). The 
expulsion of flux occurs in such a way that only part of 
the tubes leave the sample and the number of normal 
domains remain unchanged (Fig. |3jg)). The expulsion 
of the flux occurs with smaller steps compared to the 
case of magnetic field sweep up. Further decreasing the 
field leads to the formation of patterns containing singly- 
quantized vortices (Fig. [3^h)). At H = 0.023i7 c the mag- 
netic field is totally expelled and the system transits to 
the Meissner state. Notice that the square symmetry of 
the IS flux structures is always preserved, while in large 
samples the square symmetry is broken when decreasing 
the field [24] . This difference is a consequence of stronger 
confinement (surface barrier) in mesoscopic samples. 

From Fig. [5|a) it is clear that the ground state of the 
sample is not reached during both field sweep up and 
down (except for the Meissner phase) - the system is 
locked in higher energy metastable states. The ground 
state can be obtained during field cooling simulations, the 
results of which are shown by symbols in Fig. [2] As seen 
from this figure the system can nucleate in many different 
(meta) stable states when starting from different random 
initial conditions. The total number of flux quanta TV 
inside the sample is smaller during magnetic field sweep 
up (solid circles) and it is much larger when decreasing 
the field (open circles) as compared to the one in the 
ground state (stars) (see the inset of Fig. [2ja)). The 
existence of such superheated and supercooled states re- 
sults in a large hysteresis in the magnetization curve [25] . 
The magnetization of the sample, defined as the flux ex- 
pelled from the sample M = ((H) — H)/Ait, where (H) 
is averaged over the sample volume magnetic field, [20] 
is shown in Fig. [2^b). Each jump in the magnetization 
curve corresponds to a transition between different flux 
states. The difficulty of expelling flux results in a posi- 
tive magnetization in decreasing field (thin blue curve in 
Fig. [2^b)). The hysteresis disappears at very small fields, 
and at larger magnetic fields when only surface super- 
conductivity exists. Transition from the tubular state to 
the laminar-like state (4) occurs at H~QA2H C , which is 
slightly larger than the field (H c± 0.3H C ) corresponding 
to the transition between tubular and laminar states in 
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FIG. 4: (Color online) The same as Fig. [5] but for a sphere of 
radius R = 12.4£. Isosurface plots show the ground state flux 
structures in the sphere at the magnetic fields: H = 0.23H C 
(N=3) (1), H = 0.31F C (N=10) (2), H = 0.51H C (N=38) 
(3), H = 0.59H C (N=46) (4), H = 0.62H C (N=48) (5) and 
if = 0.73H C (N=55) (6). Light (dark) blue corresponds to 
superconducting (normal) region. Inset 7 shows the magnetic 
field distribution around the sphere corresponding to the state 
shown in inset 5. 



bulk samples [14]. The corresponding ground state flux 
structures are shown in Figs. [3ji-l). 

Next, we contrast the results of the cube to the case 
of a mesoscopic sphere. Figure [4] shows the free energy 
(a) and magnetization (b) vs. the magnetic field H for 
a sphere of radius R = 12. 4£. The effect of the sample 
geometry leads to the following differences. First, the 
Meissner effect becomes more pronounced leading to a 
larger first penetration field (H p = 0A7H C ) (thick red 
curve) [2TJ [26] . Due to the nonzero demagnetization fac- 
tor of the sphere the magnetic field is strongly disturbed 
near the sample boundary, resembling the magnetic field 
profile of a dipole (see inset 7 in Fig. [3]). The jump in 
the total number of flux quanta becomes larger as com- 
pared to the case of the cube (compare filled circles in 
the insets of Figs. [2|a) and[4^a)). Second, the number 
of possible meta(stable) states (hyphens in Fig. (4^a)) is 
substantially smaller in field cooling simulations, espe- 
cially at low and high fields. Only at midrange fields a 
larger number of met ast able states is found, but energy 
levels are very close to each other. The most interest- 
ing results are found for the magnetic field sweep down 
regime: the jumps in the free energy (and magnetiza- 



tion) curve (thin blue curve) are much smaller, leading 
to an almost continuous change in the number of flux 
quanta. Moreover, the system follows very closely the 
ground state, with only a few more flux quanta trapped 
inside the sample (compare open circles and stars in the 
inset of Fig. |3fa)). As a result of all this, the magnetic 
hysteresis is smaller and no paramagnetic effect is ob- 
served (Fig. |4jb)). These results allow us to conclude 
that the surface barrier for penetration of magnetic field 
in a sphere is larger compared to the one in flat samples 
(i.e., a cube), while it is strongly reduced during flux ex- 
pulsion. However, in contrast to the macroscopic sphere 
P2], there is no complete disappearance of the surface 
barrier in our mesoscopic sphere. Isosurface plots (1-6) 
in Fig. [4] shows the ground state flux structures of the 
sphere for different magnetic field values. At small fields 
only singly quantized vortices are found in the ground 
state (inset 1) and, with increasing H, giant vortices ap- 
pear (inset 2). Due to the boundary condition the flux 
tubes should be perpendicular to the surface of the sam- 
ple, leading to the curvature of flux tubes. Laminar-like 
structures first appear in the form of deformed flux tubes 
(inset 3) and radially oriented stripes (inset 4) or a ring 
of normal domain (inset 5). And, finally, we arrive at 
the surface superconductivity state with one big central 
normal domain (inset 6). 

Concluding, a remarkable variety of possible flux struc- 
tures in the IS of type-I superconductors are found in the 
mesoscopic regime. For example, single-quantized flux 
tubes (vortices) are stabilized due to the confinement im- 
posed by the sample surface. The latter also leads to the 
formation of a ring of normal domain at larger applied 
fields. All these findings are summarized into the phase 
diagram (Fig. [I]). The transition point between differ- 
ent flux structures depends on the geometry of the sam- 
ple which is related to the surface barrier for flux entry 
and exit. Regardless of the sample geometry the ground 
state of the IS at low fields consists of flux tubes and the 
transition to the laminar-like state takes place at larger 
magnetic fields. 
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